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ABSTRACT 

Transit timing variations provide a powerful tool for confirming and characterizing transiting plan- 
ets, as well as detecting non-transiting planets. We report the results an updated TTV analysis for 
1481 planet candidates (Borucki et al. 2011; Batalha et al. 2012) based on transit times measured 
during the first sixteen months of Kepler observations. We present 39 strong TTV candidates based 
on long-term trends (2.8% of suitable data sets). We present another 136 weaker TTV candidates 
(9.8% of suitable data sets) based on excess scatter of TTV measurements about a linear ephemeris. 
We anticipate that several of these planet candidates could be confirmed and perhaps characterized 
with more detailed TTV analyses using publicly available Kepler observations. For many others, 
Kepler has observed a long-term TTV trend, but an extended Kepler mission will be required to 
characterize the system via TTVs. We find that the occurrence rate of planet candidates that show 
TTVs is significantly increased (~68%) for planet candidates transiting stars with multiple transiting 
planet candidate when compared to planet candidates transiting stars with a single transiting planet 
candidate. 

Subject headings: planetary systems; planets and satellites: detection, dynamical evolution and sta- 
bility; techniques: miscellaneous 



1. INTRODUCTION 

Long-term, high-precision photometric observations 
from NASA's Kepler mission provide the basis for con- 
firming and characterizing planets via transit timing vari- 
ations (TTVs). Kepler has confirmed 10 transiting plan- 
ets via TTVs (Holman et al. 2010; Ballard et al. 2011; 
Cochran et al. 2011; Lissauer et al. 2011a). Many more 
Kepler planet candidates are expected to reveal signifi- 
cant TTVs, as the time baseline of Kepler observations 
grows (Ford et al. 2011; hereafter Fll). Fll reported 
preliminary indications of TTVs for dozens of additional 
planet candidates and predicted significant TTVs for sev- 
eral planet candidates in multiple transiting planet can- 
didate systems based on the first fourth months of Kepler 
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observations. This paper and Steffen et al. (2012b) pro- 
vide updated TTV analyses of Kepler planet candidates 
to identify planets candidates worthy of more detailed 
TTV analyses. We describe our methods in ij^l Ta- 
ble Q] provides a description of the online-only, machine- 
readable material. We provide an overview of our results 
and discuss a subset of particularly interesting systems 
in SI 

2. METHODS 

We analyze transit times observed by Kepler through 
the end of quarter six (September 20, 2010; J. Rowe et 
al. in preparation; hereafter R12). We include all active 
planet candidates from Borucki et al. (2011) and Batalha 
et al. (2012) for which there are at least three transits 
and the photometric signal-to-noise of a typical transit 
is at least three. 

For each planet candidate to be analyzed, we calculate 
the best-fit linear ephemeris (Li) to the transit times and 
uncertainties from R12. We anticipate that the distribu- 
tion of transit timing errors is likely to approximated 
by a mixture distribution with heavily weighted Gaus- 
sian component (core) and a weakly weighted distribu- 
tion with broader tails, due to the difficulties of measur- 
ing transit times (e.g., data artifacts, stellar variability, 
nearly simultaneous transit of another planet, large tran- 
sit timing variations, or other complications). In order 
to reduce the effects of a small number of poorly mea- 
sured transit times, we do not report the initial best- 
fit linear ephemeris (Li), but rather report an updated 
linear ephemeris (L2) calculated after cleaning the tran- 
sit times to exclude outliers and transit times with un- 
usually large uncertainties. Specifically, we reject any 
transit times with an absolute deviation from the linear 
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ephemeris exceeding four times the median absolute devi- 
ation of transit times from the first linear ephemeris (Li). 
We also reject any transit times for which the measure- 
ment uncertainty exceeds twice the median of the pub- 
lished transit time uncertainties. We record the median 
timing uncertainty for the remaining Ntt transit times 
(<7tt)- Then, we calculated an updated the best-fit lin- 
ear ephemeris (L2) using the cleaned transit times. We 
report the epoch (Eu n ) and period (Pn n ) of the updated 
linear ephemeris (L2) in the electronic table. 

First, we identify planet candidates with larger than 
expected scatter about the updated linear ephemeris 
(L2) in the measured transit times based on the mea- 
surement uncertainties (hereafter "excess scatter"). We 
calculate the median absolute deviation (MAD) and 
the weighted root mean square deviation (WRMS) of 
the clipped transit times relative to the updated linear 
ephemeris (Fll). We calculate the p-value for a stan- 
dard x 2 test, p x 2 j in which is the probability that we 
would measure a x 2 as large or larger than X\i n assum- 
ing transit timing errors follow a Gaussian distribution. 
We also calculate (p^^), the p- value for a x 2 -test us- 
ing x ' 2 = 7T N TT (MAD) 2 /(2a^ T ), as described in Fll. 
Since we cannot be confident that the transit time er- 
rors follow a Gaussian distribution, one should interpret 
the listed p- values with caution (particularly for low S /N 
transits) . We consider a small p- value as indicating that 
a transit timing data set is worthy of further investiga- 
tion. We identify planet candidates with ap x 2 j in < 10~ 3 
and/or p x ,i j in < 10~ 3 by setting the 1 and 2 bits of the 
TTV flag (see electronic table). For a few planet candi- 
dates, the transit times have a small number of transits, 
but each has a highly significant deviation from a linear 
ephemeris. We indicate these (MAD> 5(Ttt) with the 
3rd bit (with a value of 4) of the TTV flag. 

Second, we identify planet candidates with long-term 
trends in the deviations of the measured transit times 
from a linear ephemeris (hereafter referred to as "long- 
term TTV trends" ) . We calculate the best-fit (minimum 
X 2 ) ^-degree polynomial ephemeris 

in = E rid +nx P rid ^1 + c i" <_1 J . (!) 

where t n is the predicted transit time of the nth transit, 
F„ d and P nd are the epoch and period of the n^-degree 
model and the Cj's are the higher-order coefficients. Note 
that we renumber the transits relative to R12, so the 
epoch falls near the middle of the observed transits, re- 
ducing covariance between the epoch and period. We 
calculate Xn d f° r each rid up to 4, if there are at least 
8 transits and otherwise for rid <nTT-2. We perform 
standard F-tests using the ratio of % 2 for the n^-degree 
polynomial model and all lower-degree ephemerides. If 
Pi? quadjin, the p- value for the F-test comparing the linear 
and quadratic models is less than 2%, then we conclude 
that the difference in Xq U ad an d Xun * s greater than one 
would expect due to mere fluctuations at the 2% signif- 
icance level. Therefore, we accept the quadratic model 
and set the 4th bit (with a value of 8) of the TTV Flag 
in the electronic table . We repeat the procedure for the 
cubic and then the quartic models, accepting the higher 



degree models if the F-test p-value is less than 2% when 
comparing the higher degree model to both the linear 
model and the last accepted model. We list iVj, the de- 
gree of the highest-degree accepted model, along with 
p-values for the various tests in the electronic table. 

3. RESULTS 

We identified 39 TTV candidates with long-term 
trends using simple polynomial models (see Figures 1, 
2 and 30 Of these, we approximate the TTVs with a 
quadratic model for 14 planet candidates (KOIs 157.04, 
168.01, 227.01, 308.01, 448.02, 473.01, 784.01, 886.01, 
1081.01, 1529.01, 1573.01, 1581.01, 1599.01, 1840.01), 
a cubic model for 14 planet candidates (KOIs 456.01, 

524.01, 738.01, 806.03, 854.01, 884.02, 918.01, 1102.02, 
1145.01, 1199.01, 1270.02, 1366.02) and a quartic model 
for 11 planet candidates (KOIs 103.01, 137.02, 142.01, 

244.02, 248.01, 277.01, 377.01, 377.02, 841.01, 984.01, 
1102.01). Already, KOIs 137.02 (Kepler-18d), 157.04 
(Kepler-llf), 168.01 (Kepler-23c), 244.02 (Kepler-25b), 
248.01, 377.01 (Kepler-9b), 377.02 (Kepler-9c), 738.01 
(Kepler-29b), 806.03 (Kepler-30b), 841.01 (Kepler-27b), 
1102.01 (Kepler-24c), 1102.02 (Kepler-24b) have been 
confirmed (Cochran et al. 2011; Lissauer et al. 2011a; 
Ford et al. 2012; Steffen et al. 2012; Holman et al. 2010; 
Fabrycky et al. 2012), most occurring between this pa- 
per's original submission and publication. 

Five planet candidates have extremely significant 
TTVs, but polynomial models did not provide a signif- 
icant improvement in fit, due to the limited number of 
transits (351.02, 806.01, 1271.01 and 1474.01) or a com- 
plex signal that was not well approximated by a < 4th 
degree polynomial (872.01). While we regard these as 
relatively strong TTV candidates, more detailed analy- 
ses will be necessary before most of the planets are con- 
firmed, as one must ensure that the TTVs are dynamical 
in origin and that we are not observing eclipse timing 
variations in a triple star system (Carter et al. 2011; 
Slawson et al. 2011; Steffen et al. 2011) that has been 
diluted so the eclipse depth is consistent with a planet 
transit. Already, KOIs 806.01 (Kepler-30d) and 872.01 
(KOI-872b) have been promptly confirmed between this 
paper's original submission and publication (Fabrycky et 
al. 2012; Nesvorny et al. 2012). 

For an additional 131 planet candidates, the measured 
transit times appear to have a scatter greater than that 
expected based on the measurement uncertainties. While 
we regard these as weaker TTV candidates, we are opti- 
mistic that many may be closely-packed multiple planet 
systems for which TTVs could characterize the planet 
masses and orbital properties. For example, the elec- 
tronic table identifies several confirmed planets as having 
excess TTV scatter but without a significant long-term 
trend (KOI 72.01=Kepler-10b, KOI 137.01=Kepler- 
18c, KOI 157.03=Kepler-lle, KOI 203.01=Kepler-17b, 
KOI 244.01=Kepler-25c, KOI 250.01=Kepler-26b, KOI 
806.02=Kepler-30c, KOI 952.02=Kepler-32c). Thus, we 
recommend these planet candidates with excess TTV 
scatter for further analysis, as well as long-term moni- 
toring by Kepler using the 1-minute cadence to improve 
the number, accuracy and precision of transit time mea- 

14 Similar figures for additional planet candidates will be avail- 
able at http://www.astro.ufl.edu/~eford/data/kepler/. 
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surements. 

We note that more detailed analyses have detected 
TTVs for several planets that were not identified as hav- 
ing a long-term trend or excess scatter by our analysis 
(e.g., Kepler-llb, Kepler-19b). Thus, we anticipate that 
more detailed analyses (particularly those that use short- 
cadence observations or systems with multiple transiting 
planet candidates) may confirm and characterize addi- 
tional planets via TTVs. 

This analysis significantly improves on that of Fll, pri- 
marily due to the increased number and time span of 
transit timing observations now available. Of those with 
an obvious TTV pattern (TTV Flag of 1 in Fll Table 5), 
five of seven are identified as significant in our updated 
analysis. One exception (KOI 928.01) was not included 
in our analysis, as it is now a known triple star (Steffen 
et al. 2011b). For the other exception (KOI 528.01), we 
now rule out a long-term TTV trend, but low amplitude 
quaesiperidic TTVs are still a possibility. The remaining 
five TTV candidates have a statistically significant long- 
term trend once we include observations through Q6. Of 
the seven candidates with PF,quad,iin 

< 0.03 from Fll 

(Table 4), only two (KOI 142.01, 227.01) were identified 
as TTV candidates in our updated analysis. As noted in 
Fll, the small number of transits and small time span of 
observations was a significant limitation. 

4. DISCUSSION 

We compare the frequency of TTV candidates accord- 
ing to orbital period, transit S/N and whether the host 
star has multiple transiting planet candidates. For this 
purpose, we restrict our attention to planet candidates 
with at least 5 transits (after clipping points with un- 
usually large uncertainties or absolute deviations, as de- 
scribed in Sj2]), leaving 1386 TTV data sets, including 851 
(535) planet candidates associated with KOIs that have 
a single (multiple) transiting planet candidates. We find 
163 (18%) TTV data sets with some indication of po- 
tential TTVs (TTV Flag>l), including 89 (10.5%) and 
74 (13.8%) for stars with single and multiple transiting 
planet candidates, respectively. If we restrict our at- 
tention on systems with more robust signature of TTVs 
(TTV Flag>2), then we find 109 (7.9%) TTV data sets 
had excessive scatter including 57 (6.7%) and 52 (9.7%) 
for single and multiple candidate systems, respectively. 
Based on Monte Carlo simulations, the probability of 
these two subsets of KOIs have the same occurrence rate 
of long term TTV trends is 0.7%. 

We find 39 (2.8%) TTV data sets with a significant 
long-term trend (N d >= 2), including 19 (2.2%) and 
20 (3.7%) for stars with single and multiple transiting 
planet candidates, suggesting a ^68% greater rate of 
long-term TTV trends in systems with multiple tran- 
siting planet candidates. Based on Monte Carlo simula- 
tions, the probability of these two subsets of KOIs have 
the same occurrence rate of long term TTV trends is 
2.6%. Therefore, there is significant evidence that plan- 
ets in systems with multiple transiting planet candidates 
are significantly more likely to show TTVs than planets 
in stars with a single known transiting planet. This su- 
percedes the results of a similar analysis performed using 
just the first four months of Kepler data in Ford et al. 
(2011) that did not find a difference in the frequency of 
TTV signals between planet candidates transiting stars 



with one or multiple transiting planet candidates. 

Such comparisons provide information about the fre- 
quency of non-transiting planets with similar orbital dis- 
tances as the transiting planets and can help break the 
degeneracy between planet abundance and the distri- 
butions of mutual inclinations (Lissauer et al. 2011b, 
Tremaine & Dong 2011). Since we do not know the distri- 
bution of non-transiting planets in these systems, we can 
not clearly identify whether this difference is due to ob- 
servational bias. For planetary systems with small, but 
non-zero mutual inclinations, the probability of multi- 
ple planets transiting is increased for closely-spaced sys- 
tems that are also likely to yield large TTVs. If such 
systems are common, then one would expect the rate 
of detectable TTVs to be greater when there are mul- 
tiple transiting planets. On the other hand, the geo- 
metric probability of multiple planets transiting is dras- 
tically reduced for systems with modest mutual inclina- 
tions (> 5°; see figure 13 of Lissauer et al. 2011b). If the 
typical planetary system contained closely-spaced plan- 
ets with mutual inclinations ~ 5° or larger, then one 
would expect a similar rate of TTVs, regardless of the 
number of transiting planets, since the overwhelming ma- 
jority of multiple planet systems would have only a single 
transiting planet. Therefore, the measurement of an in- 
creased rate of TTVs for stars with multiple transiting 
planet candidates is important because it provides in- 
formation about the distribution of mutual inclinations. 
A complete analysis accounting for all relevant detec- 
tion biases is important. For example, the magnitude of 
TTVs can increase rapidly with increasing orbital eccen- 
tricity. If there were a correlation between eccentricities 
and mutual inclinations, then the increased eccentricity 
of systems with large inclinations would increase the rate 
of detectable TTVs for multiple-planet systems with only 
one transiting planet. A complete analysis accounting for 
all relevant detection biases will be quite involved and is 
beyond the scope of this paper. 

The true frequency of TTVs is likely to rise for both 
samples with more detailed analyses and as the time span 
of Kepler observations increases. Further analysis will 
be required to understand how these results are affected 
by potential complications (e.g., some multiple planet 
systems might be overlooked due to TTVs interfering 
with their detection or transit time measurement). 

We divide the TTV data sets into two equal-sized sam- 
ples (693) according to orbital period. This results in di- 
viding the sample at a period of 15. Od. We find no indica- 
tion of a difference in the frequency of planet candidates 
with some indication of potential TTVs (TTV Flag>l) 
among planet candidates with orbital period greater or 
less than 15.0d (12.1% and 11.4% for the short and long 
period sub-samples, respectively). If we restrict our at- 
tention to planet candidates with a more robust indi- 
cation of potential TTVs (TTV Flag>2), then we find 
a ^60% greater rate of TTV candidates for the long- 
period sample (6.1% and 9.7% for the short and long pe- 
riod sub-samples, respectively). Based on Monte Carlo 
simulations, the probability of these two subsets of KOIs 
have the same occurrence rate of TTV signals is <0.1%. 
However, we find no significant difference for the rate 
of long-term TTV trends among planet candidates with 
orbital period greater or less than 15. Od (i.e., 2.6% and 
3.0% for the short and long period sub-samples, respec- 
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tively). It would not be unexpected for planet candidates 
with larger orbital periods to have a greater frequency 
of detectable TTVs, since; the TTV amplitude typically 
scales with the orbital period. However, our results do 
not provide strong support for such a trend. 

We divide the TTV data sets into two equal-sized sam- 
ples according to photometric S/N of the transits. This 
results in dividing the sample at a S/N of individual tran- 
sits of 6.6. We find ~ 22% more planet candidates with 
some indication of potential TTVs (TTV Flag>l) among 
the high S/N sample and ~ 73% more planet candidates 
with a more robust indication of potential TTVs (TTV 
Flag>2). We find a <~ 37% larger rate of planet candi- 
dates with long-term TTV trends among the high S/N 
subsample. The increased frequency of planet candidates 
with excess TTV scatter or long-term TTV trends is 
likely due to the higher quality transit timing precission 
for planet candidates with high S/N per transit. 

We recommend that observers planning follow-up ob- 
servations of Kepler planet candidates check whether 
there are indications of significant TTVs. While most 
KOIs show TTVs that are small compared to the tran- 
sit duration, some systems (e.g., KOI 806) have TTVs 
so large that they can affect the scheduling of follow-up 
observations (e.g., Tingley et al. 2011; Fabrycky et al. 
2012). Based on synthetic planet populations (Lissauer 
et al. 2011b), we anticipate that ~ 15 — 26% of plan- 
ets could have TTVs as large as the transit duration. 
Of course, the higher-order ephemerides should not be 
trusted to predict future transit times accurately, as ex- 
trapolations are often unreliable. 

On sufficiently long timescales, most TTV are likely 
better modeled with periodic functions. Indeed, Steffen 
et al. (2012) have performed a complementary analysis 
using a sinusoidal model. A periodic model is likely su- 
perior for TTV signals once they have been observed 
for multiple TTV cycles. We caution TTV signatures 
can be quite complex (e.g., Veras et al. 2011), so one 
can not presume that a sinusoidal model is necessarily 
a better model until the TTV signal has been observed 
for multiple cycles. Further, polynomial models require 
fewer model parameters and do not require searching 
over many possible TTV frequencies, making polynomial 
models more sensitive to low-amplitude, long-term TTV 
signals. Thus, polynomial and sinusoidal TTV models 
are both important and complementary tools to search 
for TTV candidates. We anticipate that the sinusoidal 
model is likely more appropriate for TTV candidates that 
require a Nd > 4 degree model. 

The largest planetary TTV signatures are likely to 
arise in systems in or near mean motion resonances, 
but the associated timescales are typically multiple years 
(Agol et al. 2005; Holman & Murray 2005; Veras et al. 
2011; Veras & Ford 2011). Over timescales short com- 
pared to the TTV timescale, polynomial models can of- 
ten approximate the dominant TTV signal with a mini- 



mal number of model parameters. On the other hand, a 
periodic TTV model requires at least five model param- 
eters (and a minimum of six transits to assess a model). 
Searching over many possible TTV frequencies requires 
performing many statistical tests, reducing the sensitiv- 
ity to small amplitude signals when only a single TTV 
cycle has been observed. For planets with orbital peri- 
ods approaching a year, it will be important to analyze 
TTV signals using as few transits as possible. Borucki 
et al. (2012) obtained early TTV constraints for a planet 
in the habitable zone using just 3 transit times. For 
tightly-packed systems with multiple transiting planets, 
a full TTV analysis may be possible even with only ^6 
transits per planet (e.g., Holman et al. 2010; Fabrycky 
et al. 2012). Thus, we anticipate that polynomial TTV 
models will remain an important tool for assessing TTV 
candidates, particularly for planet candidates with long 
orbital periods, such as those approaching the habitable 
zone of solar-type stars. 

We anticipate that many of the planet candidates with 
long-term TTV trends (see Fig. 1, 2 & 3) are likely to be 
confirmed with additional analysis. Such large and co- 
herent trends are unlikely to be due to astrophysical noise 
or non-Gaussian measurement errors. In some cases, fur- 
ther observations and/or TTV analysis may lead to the 
detection of additional non-transiting planets (e.g., Bal- 
lard et al. 2011), or even small planets with large TTVs 
that were not detected by the standard Kepler planet 
search pipeline. In other cases, long-term TTV trends 
may be due to perturbations from stellar or brown dwarf 
companions rather than additional planets. It is also 
possible that a transiting planet candidate with TTVs is 
actually an eclipsing binary that is significantly diluted, 
so the transit depth is consistent with planetary transit, 
and a member of a triple star system (e.g., KOI 928.01; 
Steffen et al. 2011). It is less clear what fraction of the 
planet candidates with excess TTV scatter have dynam- 
ically induced TTVs, as opposed to an apparent TTV 
signal due to astrophysical noise (e.g., Desert et al. 2011) 
or non-Gaussian measurement errors. Nevertheless, we 
recommend that these candidates be subjected to more 
detailed TTV analysis, as some of these are planets in 
closely-packed systems with complex TTVs (e.g., Kepler- 
11; Lissauer et al. 2011a). Such systems are particularly 
favorable for measuring planet masses and bulk density, 
so as to constrain the planet's composition and nature 
(i.e., distinguish between rocky versus gaseous planets). 



Funding for this mission is provided by NASA's Science 
Mission Directorate. We thank the entire Kepler team 
for the many years of work that is proving so successful. 
E.B.F acknowledges support by the National Aeronautics 
and Space Administration under grant NNX08AR04G is- 
sued through the Kepler Participating Scientist Program. 
Facilities: Kepler. 



REFERENCES 



Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 
359, 567 

Ballard, S., Fabrycky, D., Fressin, F., et al. 2011, ApJ, 743, 200 

Batalha, N. M., ct al. 2012, to be submitted to ApJS. 

Borucki, W. J., Koch, D. G., Basri, G., et al. 2011a, ApJ, 736, 19 



Borucki, W. J., Koch, D. G., Batalha, N., et al. 2012, ApJ, 745, 
120 

Carter, J. A., Fabrycky, D. C, Ragozzine, D., et al. 2011, Science, 
331, 562 

Cochran, W. D., Fabrycky, D. C, Torres, G., ct al. 2011, ApJS, 
197, 7 



TTV Candidates from Kepler 



5 



Desert, J.-M., Charbonneau, D., Demory, B.-O., et al. 2011, 
ApJS, 197, 14 

Fabrycky, D. C, Ford, E. B., Steffen, J. H., et al. 2012, ApJ, 750, 
114 

Ford, E. B., Rowe, J. F., Fabrycky, D. C, et al. 2011, ApJS, 197, 
2 

Ford, E. B., Fabrycky, D. C, Steffen, J. H., et al. 2012, ApJ, 750, 
113 

Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288. 
Holman, M. J., Fabrycky, D. C, Ragozzine, D., et al. 2010, 
Science, 330, 51 

Lissauer, J. J., Fabrycky, D. C, Ford, E. B., et al. 2011a, Nature, 
470, 53 



Lissauer, J. J., Ragozzine, D., Fabrycky, D. C, et al. 2011b, 

ApJS, 197, 8 
Nesvorny, D., et al. 2012, Science, 336, 1133. 
Slawson, R. W., Prsa, A., Welsh, W. F., et al. 2011, AJ, 142, 160 
Steffen, J. H., Quinn, S. N., Borucki, W. J., et al. 2011, MNRAS, 

417, L31 

Steffen, J. H., Fabrycky, D. C, Ford, E. B., et al. 2012a, MNRAS, 
421, 2342 

Steffen, J. H., et al. 2012b, ApJ, in press. larX!v:1201. 18731 
Tingley, B., Palle, E., Parviaine n, H., et al. 2011 , A&A, 536, L9 
Tremaine, S., & Dong, S. 2011, arXiv:1106.5403 
Veras, D., Ford, E. B., & Payne, M. J. 2011, ApJ, 727, 74 
Veras, D., & Ford, E. B. 2011, MNRAS, L369 



6 Ford et al. 



TABLE 1 

Kepler TTV Metrics Catalog Format 



Column 


Format 


Description 




1 


KOI 


F7.2 


Kepler Object of Interest Identifier 


2 


S/N 


F7.1 


Photometric signal-to-noise for typical transit 


3 


nTT 


15 


Number of transits used for analysis 


4 




F5.1 


Estimated transit duration (Batalha ct al. 2012) 


5 




12 


Degree of polynomial selected for ephemeris 


6 


r I If I IT T T711 

1 1 V b lag 


14 


1 1 V t lag 


I 


rp 

-c-lin 


r LZ.O 


Epoch of best-fit linear ephemeris 


o 


°\E lin 


r y.u 


Uncertainty in E\^ n 


9 


Pi- 
1 lin 


F11.6 


Period of best-fit linear ephemeris 


10 


ffp,. 


F9.6 
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Median transit time uncertainty in minutes 
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Median absolute deviation of transit times from a linear ephemeris in minutes 


23 


WRMS 


F6.1 


Weighted root mean square deviation of transit times from a linear ephemeris in minutes 
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p-value for a x 2 -test using x' 2 relative to quadratic ephemeris 
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p-value for a x 2 -test relative to cubic ephemeris 
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p-value for a x 2 -test using x' 2 relative to cubic ephemeris 
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p-value for a x 2 -test relative to quartic ephemeris 
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p-valuc for a x 2 -test using x' 2 relative to quartic ephemeris 
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p-valuc for an P-test comparing linear and quadratic ephemerides* 
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p-value for an P-test comparing linear and cubic ephemerides* 
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p-value for an P-test comparing quadratic and cubic ephemerides* 
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p-value for an P-test comparing linear and quartic ephemerides* 
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p-value for an P-test comparing quadratic and quartic ephemerides* 
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p-value for an P-test comparing cubic and quartic ephemerides* 



Note. — All ephemerides and TTV statistics are calculated from Kepler observations up to and including quarter six (Rowe et al. 2012). 
a Sum of following flags: 1 if p %2 < 10~ 3 , 2 if p x , 2 < 10~ 3 , 4 if MAD > 5<r TT , 8 if PF.lin.quad < 0.02 AND nTT> 4, 16 if PF.lin.cubc < 0.02 
AND nTT> 5 AND IF(8-bit flag is set, PF, q uad,<mbc < 0.02, TRUE), and 32 if PF,ii„, q uar < 0.02 AND nTT> 8 AND IF(16-bit flag is set, 
PF.cubc.quar < 0.02, IF(8-bit flag is set, p F , quad , q uar < 0.02, TRUE)). 
b Zero if N d <2 
° Zero if N d <3 
d Zero if N d <4 

° For description of x' 2 statistic and tests, see § 3.1 of Ford et al. 2011 
* Value of -1.0 if insufficient number of transits 
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Fig. 1. — Transit times for planet candidates with long-term TTV trends. Kepler-9c (panel c), Kepler-30b (panel d), Kepler-30d (panel 

g) , Kepler-9b (panel i) and Kepler-24b (panel 1) have already been confirmed via TTVs (Holman et al. 2010; Fabrycky et al. 2012; Ford et 
al. 2012). Dotted lines show the best-fit polynomial model given in Tablc[l] For KOIs 351.01 (panel f), 806.01 (panel g) and 1271.01 (panel 

h) there is no polynomial model, since a polynomial model (of degree 4 or less) did not significantly improve upon the linear ephemeris. 
Similar figures for additional planet candidates will be available at http://www.astro.ufl.edu/~eford/data/kepler/. 
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Fig. 2. — Transit times for planet candidates with long-term TTV trends. KOI 872b (panel e), Kepler-24c (panel f), Kepler-23c (panel 
i) and Kepler-29b (panel o) have already been confirmed via TTVs (Nesvorny et al. 2012; Ford et al. 2012; Fabrycky et al. 2012). Dotted 
lines show the best-fit polynomial model given in Table [l] For KOI 872.01 (panel e) there is no polynomial model, since a polynomial 
model (of degree 4 or less) did not significantly improve upon the linear ephemeris. 
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Fig. 3. — Transit times for planet candidates with long-term TTV trends. Kepler-27b (panel b), Kepler-llf (panel f), Kepler-18d (panel 
m) and Kepler-25b (panel n) have already been confirmed via TTVs (Steffen et al. 2012; Lissauer et al. 2011a; Cochran et al. 2011). Dotted 
lines show the best-fit polynomial model given in Table [T] For KOI 1474.01 (panel d) there is no polynomial model, since a polynomial 
model (of degree 4 or less) did not significantly improve upon the linear ephemeris. 



